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Abstract 

The active cloak comprises a discrete set of multipole sources that destructively 
interfere with an incident time harmonic scalar wave to produce zero total field over 
a finite spatial region. For a given number of sources and their positions in two di- 
mensions it is shown that the multipole amplitudes can be expressed as infinite sums 
of the coefficients of the incident wave decomposed into regular Bessel functions. The 
field generated by the active sources vanishes in the infinite region exterior to a set of 
circles defined by the relative positions of the sources. The results provide a direct 
solution to the inverse problem of determining the source amplitudes. They also define 
a broad class of non-radiating discrete sources. 

1 Introduction 

Cloaking is intended to make an object undetectable to incident waves. The approaches 
proposed consist mainly of two quite distinct types of cloaking, namely passive and active. 
Passive cloaking requires devising a metamaterial that can steer the wave energy around the 
object regardless of the incident wave. Our interest here is with active cloaking, specifically 
in situations where the active sources lie in the exterior of the region containing the cloaked 
object. We call this configuration active exterior cloaking in keeping with prior terminology 

m 

Despite the dominant interest in passive cloaking devices, active exterior cloaking has 
been investigated quite extensively [21 [U [3l HJ |5]. Miller [2j proposed creating a cloaking 
region by measuring particle motion near the surface of the cloaking zone while simulta- 
neously exciting appropriate surface sources where each source amplitude depends on the 
measurements at all sensing points. As an active cloaking method this approach is limited 
because it does not provide a unique relationship between the incident field on the one hand, 
and the source amplitudes on the other. A solution to this problem was provided by Vasquez 
et al. [H E] in the context of active exterior cloaking for the 2D Helmholtz equation. They 
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used Green's formula and addition theorems for Bessel functions to formulate an integral 
equation, which was converted to a linear system of equations for the unknown amplitudes. 
Crucially, the integral equation provides the source amplitudes as linear functions of the 
incident wave field. Vasquez et al. also showed, by construction, that active cloaking can be 
realized using as few as three active sources in 2D. A more explicit form of the linear relation 
for the source amplitudes as a function of the incident field was developed in f5] . Multipolar 
sources were used to reproduce Miller's cloak and numerical results were compared with 
S VD solutions of the linearized system jU [3] . The approach of [U [3] was generalized in [4J to 
handle the 3D Helmholtz equation, seeking non resonant frequencies of the cloaked object. 
Further analysis and extension of the methods to the quasistatic regime relevant to Laplace's 
equation can be found in P, El E]. The active source method of Vasquez et al. has also been 
adapted to create illusion effects so that an object outside the cloaking region can be made 
to appear like another object [9j. 

In this paper we demonstrate that the integral representation of Vasquez et al. [5j for the 
source amplitudes can be reduced to closed- form explicit formulas. This obviates the need 
to reduce the integral equation of Vasquez et al. [H [3] to a system of linear equations, which 
must then be solved numerically, or to evaluate line integrals, as proposed in f5]. We provide 
analytical expressions for the source amplitude coefficients for general incidence as well as 
plane wave incidence. The expressions involve no more than sums of cylinder functions which 
can be truncated to achieve any desired accuracy. We also prove that the field generated by 
the active sources vanishes in the infinite region exterior to a set of circles defined by the 
relative positions of the sources. The active source field, by construction, cancels the incident 
field in the cloaked region, which is defined by the region interior to the same circular areas. 
The analytical results are verified by calculation of the farfield and the nearfield amplitudes, 
which are shown to vanish when the summation is accurately evaluated. 

The non-radiating nature of the active field has relevance to the inverse source problem 
pH] . Although for this problem, some uniqueness results are available for restricted forms 
of sources, e.g. "minimum energy sources" [JJJ, in general the solution to this problem is 
know to be non-unique [l2] . Here we develop the solution of the active cloaking problem as 
a new family of non-radiating sources, with the property that they cancel a given incident 
field over a finite region. 

We begin in ^ with a statement of the problem, a review of the governing equations, 
and a summary of the main results, given in eqs. (Ej) and ([7]). The basic integral relation of 
Vasquez et al. is derived in ^ from which the main results are shown to follow. Some 
example applications of the new formulas are presented in ^ Some implications of the 
general results are discussed in ^ and conclusions are given in ^ 

2 The problem and its solution 
2.1 Problem overview 

The active cloaking devices considered here operate in two di mensi ons, and consist of arrays 
of point multipole sources located at positions G M^, m = 1, M, see Figure [H The active 
sources lie in the exterior region with respect to the cloaked region C and this type of cloaking 
may therefore be called "active exterior cloaking" [Ij. Objects are undetectable in the cloaked 
region by virtue of the destructive interference of the sources and the incident field with the 
result that the total wave amplitude vanishes in the cloaked region C. The advantages of this 
type of cloaking device are: (i) the cloaked region is not completely surrounded by a single 
cloaking device; (ii) only a small number of active sources are needed; (iii) the procedure 



Active exterior cloaking 



July 23, 2012 



3 




Figure 1: Insonification of the actively cloaked region C generated by M active point mul- 
tipole sources at x^, m = 1, M. The region R is defined as the interior of the union of the 
dashed circular arcs, that is, the combined area comprising C and the M circular domains. 
The incident field in this case is a plane wave with wave vector k in the direction ip. 

works for broadband input sources. A disadvantage of the active cloaking approach is that 
the fields near the ideal sources may become uncontrollably large. Realistically these would 
be replaced by regions of finite extent and thus their magnitude is reduced. A further 
disadvantage of the method is that the incident field must be known. However we note that 
with the approach proposed in this paper, the new expressions require only the expansion of 
the incident field into entire cylindrical waves, as compared with the line integrals derived 
in [4j which require knowledge of the incident field and its normal derivative. 

The shaded region in Figure [1] denotes the cloaked zone C generated by M active point 
multipole sources. The boundary of C is the closed concave union of the circular arcs m = 

1, M, {a^, (f)^^\ 02^^} associated with the source at x^. In the general case {a^, (/>i^\ 02^^} 
are distinct for different values of m. Note that the wave incidence shown in Figure [H is 
a plane wave although the solution derived below is for arbitrary incidence. The inverse 
problem to be solved is to find the amplitudes of the active sources as a function of the 
incident wave, and to prove that the cloaked region is indeed the closed region C. 

2.2 Basic equations 

We assume time harmonic dependence e~^^^ which is omitted hereafter and consider the 
scalar Helmholtz equation in two dimensions. Thus the method proposed here is applicable 
to any physical situation described as such. For ease of discussion however let us consider 
the case of acoustics, so that the governing equation for the (time harmonic) pressure u{x.) 
is 

V^u + k\ = s, (1) 

where k = uj/c is the wavenumber, c the acoustic speed, and the term s represents sources. 
For a given incident wave we assume there is an additional field resulting from the active 
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sources which exactly cancels the incident wave in some bounded regio n C. This additional 
wave field is caused by the M multipole sources located at x^, m = 1,M. The assumed 
form of the total field the incident wave and the active source field Ud are, respectively 

u = Ui + (2a) 

CO 

Ui= A^U+i^), (2b) 

M oo 

^d='^ Yl - Xm), (2c) 

m=l n=—oo 

where the wave functions U^{'x) and V^{'x) are defined by 

C/„±(x) = J„(fc|x|)e±--^^ i;±(x) = i/W(fc|x|)e±--s^. (3) 

Here argx G [0, 27r) and arg(— x) = argxivr G [0,27r). Define the derivative functions 
U^^r) as 

Ufi^) = j;(Me^^"^^^^. (4) 

In the following we write Uq and Vq, with obvious meaning. Note that the functions U^{'x) 
and V^{'x) possess the properties 

t^n^(-x) = (-I)"t/.^(X), K^(-X) = (-I)X^(X). (5) 

The active source field in (l2cl) is of the same form as considered by Vasquez et al. [3l 
eq. (5)]. The three dimensional analog is given in [4, eq. (40)]. The coefficients A^^ which 
define the incident field, include as a special case plane wave incidence in the direction 

The active cloaking problem is now to find (i) the coeflScients bm,n such that the total 
field u vanishes inside some compact region (7, and (ii) to define the region C. 



2.3 Summary of main results 

The principal results can be summarized in two theorems. The first provides necessary and 
sufficient conditions on the source amplitudes b^^n in order to ensure cloaking in the region 
C and a non-radiating source field u^- The second provides the explicit expressions for the 
active source amplitudes. 

Theorem 1. Necessary and sufficient conditions on the active source coefficients b^^i in 
order to ensure zero total field (ui + Ud = 0) inside C and no radiated field (ud ^ in the 
far field) are 

M oo I U^-li'^m) 0, 

VnGZ: ^-'^ \ (6) 

These identities provide a useful means to quantify error in active cloaking as will be 
seen later on. We now state the explicit form for the source amplitudes, together with the 
shape of the cloaked region C and the region in which the source field vanishes. 
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Theorem 2. Given M active sources located at x^, m = 1,M, the required active source 
amplitude coefficients for the general incidence (l2bl) are 



bm,i = ^^M^^ "'^^'"6 (7a) 



n=—oo 



p 



kcL f 1) 

bm,in = ^^ ^ Un+pi^m) ^ _^ [jp{ka^) J[{ka^) - Jp{kam)Jiikar, 

p=—oo ^ 

X [e-^C+f)'^^'"' _e-^(^+f)'^i'"']. (7b) 

This ensures cloaking (zero total field) in the region C which is the closed and bounded 
domain formed by taking its boundary as the closed concave union of the circular arcs defined 

by {a^, (/>i^\ 02^^} tt^rf denoted as dCm, see FigureUi These coefficients also ensure that the 
radiated field from Ud is identically zero in the region exterior to all of the circles centered at 
the source points: 

Ud{^) = for X G R^/R, R = C[j{^ : |x - x^| < a^, m = TjT}. (8) 
This is the exterior to the union of the dashed circular arcs in Figure Ui 



An alternative and more concise formulation of eq. (1751) is obtained using the notation 
of eq. © with af ) = a^e{(j)^), {t = l,2), 

bm,in = jka^ E Un+pM)-^ [[/" (a)f/r'(a) - C/-'(a)[/r(a)] (9) 

p=—oo 

where e((/)^) is a unit vector subtended at angle 0^, as illustrated in Figure HI 

An important case for which the summation in (17al) can be simplified is plane wave 
incidence. Assuming the incident field is a unit amplitude plane wave in direction Ui = 
defined by = i^e~^^^, results in 

kam ^ i^e'P"^ r , , . 4"^ plane wave . . 

&m,/ = ^^(xj^ 2^ [U^{e)U, (a)-[/^ (a)[/, (a)J ^^^^ , i^HHmce. ^^^^ 



p=—oo 



The form of the coeflBcients h^^i is discussed further below. Note that the term in f ubl ). 
(I9j) and in f lTOl ). corresponding to p + / = is zero, which follows from THopital's rule, or 
otherwise. 

Theorems H] and 12 are proved in the next section. 



3 Proofs of Theorems 1 and 2 

3.1 Theorem 1: Necessary and sufficient conditions on the source 
amphtudes 

We first prove the constraints on the source coefficients 6^^^ given by Theorem HJ and at 
the same time show that they may be interpreted in terms of the near- and far-field of the 
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active sources. To this end, we express in two different forms using the generahzed Graf 
addition theorem [13^ eq. (9.1.79)], 



Let us first consider the radiated field, assuming that Ud does not radiate energy into the 
far field. The first of f fTTI) . for |x| > |y|, allows us to rewrite as a sum of multipoles at 
the origin: 

oo 

F„]/ + (x) for |x| > max(|x^| +a^), (12) 

n=—oo 

where 

M oo 

= E E bm,iU-_,{^^). (13) 

m=l l=—oo 

Define the farfield amplitude function /(^), = argx, such that 

«d(x) = m^^^ + 0((fc|x|)-3/2), |x|^oo. (14) 
The farfield amplitude function follows from the asymptotic form of the Hankel functions as 

oo ^ 1/9 

/(^)=5]^e-^ /n=(^) (15) 

n=—oo 

A measure of the nondimensional total power radiated by the sources is given by the non- 
negative far-field fiux parameter 

ar= / de\f{e)\' = 4 \Fn\'- (16) 

n=-oo 

Since Ud does not radiate energy into the far field, the active sources must vanish, so that 
Fn = Vn. Imposing this in (fT3l) ensures the necessity of (I6j)i. The sufficiency of ([6j)i is seen 
immediately by substituting (|6j)i into f fT3l ) and f fT2l ) which gives Ud = 0. 

Now let us consider the near-field inside the cloaked region C where we assume that the 
cloaked region contains the origin and the total field is zero inside (7, i.e. Ui + Ud = 0. Using 
the second identity in f fTTI ) , the active source field can be expressed in a form that is valid 
in the neighborhood of the origin (assuming |x^| > Vm), 



u 



= ^ EnU^{x.) for |x| < min(|x^| - a^), (17) 

where 



d 

n=—oo 



M oo 



= E E bm,lV-_i{^m). (18) 



m=l l=—c)G 
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The total field vanishing in some neighbourhood of the origin thus implies that + 
vanishes for every value of n. This gives rise to the necessary condition (|6j)2. Sufficiency is 
once again immediate by assuming the form ([6]) 2 and back-substituting into the forms of Ud 
and Ui above. 

Further implications of this result are explored after we complete the proof of Theorem 

El 

3.2 Theorem 2: Explicit forms for the active source ampHtudes 

The Green's function ^'(x, x') is defined as the solution of ([I]) for source s = (5(x — x^), i.e. 
5'(x,x') = — |Vo(x — x'). Consider a region D such as that depicted in Figure O chosen 
so that it does not contain any sources. We will determine the explicit form for the active 
source amplitudes together with the form of D that ensures cloaking. The latter, already 
introduced as (7, is the region depicted in Figure [H 



X2 • 




Xl 



• X4 



Figure 2: A configuration of M = 4 sources, and a region D in which the integral identity 
(linD holds. 

By assumption, both Ui and Ud satisfy the homogeneous Helmholtz equation in D (eq. 
([I]) with 5 = Vx G i^), and therefore 

/ d5(y)[^(y)a^^(y,x) -^(y,x)a^^(y)] = ^(x), v = {ui.Ud}, x G L^. (19) 

JdD 

where dD is the boundary of D depicted in Figure [2] as the union of the arcs dDj^^ m = 1, M 
and it is traversed counter-clockwise. We wish to determine the cloaked region C C D which 
is defined by its property that the total field Ui + vanishes inside (7, so that 

Ud{^) = -u,{^) = [ dS{y)[u,{y)dnVo{y-^)-Vo{y-^)dnU,{y)], x G C. (20) 

4 JdC 



Given that the boundary of C is split up, as for D into segments dCm^'rn = 1,M, we can 
use (fTTI) i . in order to write, for some xq 



00 

Vo{y - x) = yo(x - xo - (y - xo)) = ^ V;+(x - xo)t/-(y - xo) (21) 

n=—oo 
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which holds for |x — Xo| > |y — Xo|. Do this for each of the contours choosing Xq = x^ on 
each dCjn^ , so that 



. M oo 

m=l n=—oo 



\4+(x - x^) / dS^ {ui{y)dnU^ (y - x^) - t/„ (y - x^)(9„«i(y)) 

(22) 

where we require |x — x^| > |y — x^| on each contour dCm (recall that the integral is 
being considered for x G C). The minus sign in ( l22l) arises since upon expanding about the 
point x^, the counter-clockwise orientation with respect to the centre x^ is opposite to the 
counter-clockwise traversal of dC with respect to some origin inside C. Note that for this to 
hold simultaneously for all m the contours dCm must be circular arcs as depicted in Figure 
[31 Therefore we have proved that C is the region with boundary as the closed concave union 

of the circular arcs defined by {a^, (/)^^\ 02^^} and denoted as dC^nj see Figured! Finally 
using the form for Ud given in f l2cl) . we find that 



dSr, 



dCrr 



Ui{y)dnU^ (y - x^) - {y - x^)a^^^(y) 



(23) 



This agrees with ^ Eq. (8)] apart from a factor i/4 missing there. Equation ( !23ll provides 




Figure 3: The integration curve dC split into M = 3 portions dCm appropriate for the 
integral representation ( l22l) of the active source field. The cloaked central (black) region, is 
bounded by dCm^ m = 1, 2, 3. 

a direct method for calculating the multipole source amplitudes, as has been demonstrated 
numerically for different source configurations [5j. The result is not optimal, however, as it 
requires evaluation of a line integral, which can be computationally time consuming. 

The explicit formula for the source amplitudes follows from eq. ( I23l) by introducing the 
forms for the functions U~ as follows 

bm,i = --kam / d(t)e-'^'^[ui{y)fi{kam) - Ji{kam)k-^dnUi{y)]. (24) 

We see that the cloaked region C is indeed the subdomain of D in which Graf's theorem can 
be simultaneously invoked for all of the M active sources. 
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Consider plane wave incidence in the direction of the unit vector e('0), Ui = t^^(x) where 

li^(x) = e^^^W-^ (An = Te-^^^)). (25) 

Then ( l24l ) becomes, with = fc^^, 

^(^) 

= ^«^(x^)e-^^^[j/(a^)G'(a^) - MaJG' {aj] , (26) 
where the function G is defined as 

G{a)= / d(/)e^(^^^^^-^^) = V / d(/)e-^(^+^)^. (27) 

The identity e^^^^^^ = X]^-cx) "^^(^)^^^^ been used in simphfying the form of G{a). 
Performing the integration in f l27l) . we arrive at an exphcit expression for the amphtude 
coefficients 



p=—oo ^ 

(28) 



Now consider the incident field 

27r" 



/ d7/;^^(x)e^^^ = [/ + (x) (Ap = 5.^). (29) 
Jo 



It follows from integration of f l28l ) that the general form of the amplitude coefficients for the 
general incidence f l2bl ) is given by f lTbl ). 

Finally, we turn to the question of where the active source field vanishes, noting that the 
integral ( l20l) vanishes identically for field positions outside C [14j 

[ dS{y)[u,{y)dnVo{y-^)-Vo{y-^)dnU,{y)]=0, x G MVC. (30) 

4 JdC 

How does this relate to the source field Ud{'x.)7 In the course of the derivation of the coeffi- 
cients bjn,n the field Ud{'x.) was expressed in the form (l22l) for x G C The latter restriction 
on X can be removed since it is clear that eq. f l22l ) defines Ud{'x.) for all x. This is evident 
from the definition f l2cl ) and from the identity f l23l ) for 6rn,n- Equation f l30l ) therefore implies 
that Ud{'x.) vanishes at all positions outside the cloaked region for which the representation 
(I22l) holds, i.e. {x C : |x — x^| > |y — x^|, y G dCm^ m = 1, M}. This is precisely the 
region R defined in ([8]), equal to, for instance, the exterior to the colored regions in Figure 

m 

This completes the proof of Theorem El 
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Figure 4: Plane wave insonification of the cloaking region C generated by M = 8 active 
sources. 

4 Numerical examples 
4.1 Active source configuration 

We illustrate the results for plane wave incidence on configurations of the type shown in Fig. 
m The M sources are symmetrically located on a circle, with 



a^n = a, |x^| =6, = {m — l)9o m = 1, M, where Oq = 27r/M, (31) 

and by necessity, a > 6sin The circular arcs, which all have the same angular extent, are 
then defined by 



(m) , n ^ • -1 • ^\ ^ 



m = l,M. (32) 



We take a = 6sin-^ in all examples considered. Note that the cloaked region C can be 
formed by a minimum of 3 sources. A configuration with M = 8 sources is shown in Figure 
m All calculations were performed for plane wave incidence on configurations of the type 
shown in Figure [4] with varying numbers of sources, M > 3. 



4.2 Near and farfield amplitudes 

The efficiency of the cloaked region is assessed by examining the farfield and nearfield as 
functions of various parameters. If all terms in the infinite sums in eqs. f fT3l ) and f fTSl ) are 
available then the farfield is identically zero and the nearfield exactly cancels the incident 
wave, by Theorems [1] and [2l We therefore consider truncated versions of the infinite sums 
so that the farfield and nearfield coefficients, and of eqs. ( |T3l ) and ( |T8l ) respectively, 
are approximated as 



EEV, X "-' t"" v„,ez^ (33) 
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(b) Farfield, M=3, N=15 
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Figure 5: The farfield radiation amplitudes |Fi^^^^| of eq. f l32l ) for different orders of Bessel 
functions n = —10, 10, wavenumbers A; = 1,5 and truncation values = 10 and 15. The 
configuration is M = 3 multipole sources located at the distance 6=1 from the origin, with 
angle of incidence '0 = 7°. 
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(a) Farfield, M=3, k=1 



10 



10"^ 



10" 



ipapp. 



10 



10 



10"' 









'V 








N=5 

N=10 

N=15 





-10 -5 5 

n 

(b) Farfield, M=8, k=1 
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Figure 6: Dependence of the farfield amplitudes on the order n of Bessel functions 

for different values of in f l33l) (A^ = 5, 10, 15) and for different numbers of active sources: 
(a) M = 3, and (b) M = 8. The incident wavenumber is A; = 1. 



In the limit of ^ oc exact cloaking is achieved. Restricting the summation to finite 
values of A^ is equivalent to limiting the order of the active multipole sources. The behavior 
of the approximate coefficients F^^^^^ and En^^^ has implications on the accuracy of the 
cloak regardless of the type of object to be cloaked. Thus, the farfield coefficients determine 
the radiated field everywhere outside the cloak, and must necessarily be small regardless of 
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whether or not an object is being cloaked. Similarly, the total field in the cloaked region 
C must be small in order to achieve cloaking. The two conditions correspond to Fn^^^ and 
E^i^^^ + An having small values. The examples in this subsection examine the sensitivity of 
these quantities. The sources are located at 6 = 1 with plane waves incident at ^0 = 17°. 

The farfield amplitude coefficients \Fn^^\-, n = —10, 10 are depicted in Figures [5] and [6] for 
different values of the wavenumber fc, the number of sources M, and the number of terms in 
summation f l33l ). A^. It is clear from these two ffgures that the error in the farffeld coefficients 
decreases (i) as increases, (ii) as M increases, and (iii) as k decreases. The convergence 
is particularly fast as a function of A^. For instance, at A; = 1 the farfield coefficients are 
uniformly less than 10~^ for all M > 3 if > 5. Much smaller values (10~^^ or less) for 

\Fn^^\ are easily achieved for moderate values of A^, e.g. A^ = 10. 




-5 5 

n 

(a) Nearfield, M=6, k=5 
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(b) Nearfield, M=8, k=5 
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Figure 7: Variation of the nearfield amplitude coeflBcients \An + En^^\ for different values 
of the truncation size A^ in eq. ( l33l) . generated by M = 6 active sources in (a) and M = 8 
sources in (b). In all cases k = 5. 

The nearffeld amplitude coefficients \An + En^^^^l are shown in Figures [7] to [91 In contrast 
with the farfield case relatively large values of the truncation size A^ are required to obtain 
small nearfield coefficients. Figure [7] shows that A^ on the order of 100 or more is required 
to achieve accuracy comparable to the farfield coefficients. However, unlike the farfield 
amplitudes, it is found that the nearfield coefficients generally increase in magnitude with 
|n|, the order of the Bessel functions. The relatively large values of \An + En^^^^l and their 
increase with the order |n| does not necessarily mean that the total field in the nearfield 



is divergent. For instance, the top curve in Figure [71(a) indicates \Aio + E 



(app) I 
10 



=0(102), 



but this value multiplies Jio{kr)^ and, for instance, \Jio{kr)\ < 2 x 10 ^ within C. In 
other words, the increasing values of + En^^'^l with n can be balanced by the fact that 
Jn{kr) = ^(y)"^ + ... for smaU kr. 

FigurelS] shows the dependence of the nearfield coefficients on the number of sources. 
The case of the minimum number of sources, M = 3, appears to be strikingly different from 
others (M > 4). As Figure [8] indicates, adding one more source and taking M = 4 reduces 
the error from 10^ to 10~^^ for k = l,n = ±5 and from 10~^ to 10~^^ for k = = ±5. 
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Generally, as with the farfield coefficients, increasing the number of sources improves the 
accuracy of the nearfield amplitudes \An + En^^^'l. 

Finally, Figure [9] shows the nearfield dependence on the wavenumber, A; = 1,5. The 
accuracy actually improves with increasing fc, unlike the farfield case. However, it should 
be borne in mind that the nearfield coefficients multiply the terms J^(/cr), which increase in 
magnitude with k for fixed r. 



(a) M=4 




(b) M= 



\ k=1 

^\ k=2 

\ \\ k=3 


/ / 

/ 

^ t/ 


\>\ \\ k=5 / 





10 



(d) M=10 
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Figure 9: Variation of the nearfield amplitude coefficients with number of active sources 
[M = 4, 6, 8, 10) and with wavenumber (A; = 1, 5). In all cases = 130. 



The numerical results in Figures [5] through [9] show that greater accuracy is achieved 
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using more sources, which is not unexpected. For the case of M = 3, the minimum number 
required, the nearfield coefficients could be large enough to significantly diminish the cloaking 
effect. This suggests taking M = 4 might be preferable. 




-1.5 \^ -1.5 r 

-2 -2 -2 -2 

Figure 10: Absolute value of total pressure field with 4 active sources, 6 = 1, angle of 
incidence = 17°, wave number = 2, and = 60. Values above 2 in magnitude are 
clipped to make the plots visible. 




Figure 11: Real part of total pressure field with 4 active sources, b = 1^ ip = 17°, wave 
number k 10, and 60. 
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4.3 Total field 



The total field for unit amplitude plane wave incidence on configurations of active sources 
of the type defined in §4.11 is illustrated through several examples. In all cases 6=1 and 
= 17°. Figure [TOl shows the absolute value of the field for four active sources: the subplots 
provide different perspectives, indicating that the field is indeed essentially zero in the cloaked 
region C, and that the radiated field Ud is zero outside the region R. The major variation in 
the source field is within the circular regions centered on the active sources. It is found that 
the field in these regions can take very large values, and therefore, for the sake of visibility 
we truncate the plot at an arbitrary value (here = 2). Note also that the cloaked region spills 
over slightly into the circular regions. This effect is perhaps easier to see in the subsequent 
examples. 

Figure [H] considers the same M = 4 configuration of active sources at a higher frequency 
k = 10. The plots in this case show the real part of the total field, clearly illustrating the 
plane wave field in the exterior of R. The subplot on the right clearly shows that the cloaked 
region is somewhat larger than (7, extending partly into the circular regions. The number 
of modes used in Figure [H] (A^ = 60) is more than adequate to ensure convergence and 
accurate cloaking. It is more instructive to consider the effect of fewer modes, as in Figures 
[l2]to[l4l In Figure [T2l the number of modes used is on the order of the frequency, and good 
accuracy is still observed. Notice the smaller footprints of the active sources, as compared 
with Figure [HI indicating that the higher modes "fill out" the regions where is highly 
variable. Only = 5 modes are used in Figure [131 can see the deterioration of the 

cloaking effect expected with an inadequate number of multipoles. The plane wave is clearly 
evident inside the cloaked region (7, as is some scattering effects in the "shadow" zone, it 
is interesting to note that the active source footprints are reduced in size as compared with 
Figure [121 Finally, in Figure [HI we consider the effect of a larger number of active sources 
combined with a small number of modes. Comparison of Figures [13] and [14] indicates the 
tendency observed from the results of §4.21 that more active sources improves the cloaking 
effect. This is also to be expected from the discussion below in ^ which shows that for large 
numbers of sources only the lowest order multipoles play a significant role. 




Figure 12: Real part of total pressure field with 4 active sources, ip = 17°, k = 10. The 
number of modes used in the truncated sum is here limited by = 10. 
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Figure 14: The same as in Figure [T3l except now M = 7. 



4.4 Scattering examples 

Finally, we illustrate the effect of active exterior cloaking on plane wave scattering from rigid 
and soft cylinders (Neumann and Dirichlet boundary conditions, respectively). In each case 
the cylinder is circular of radius ao = 1 centered at the origin, five active sources with 6 = 4 
are used, the frequency is = 5, and the incident wave strikes at angle = 17°. 

Figures [15] and [16] compare the response from a rigid cylinder with the active cloaking 
turned on and turned off. The absolute value is shown in Figure [T5] while Figure [T6] considers 
only the real part of the complex field, which clearly indicates the plane wave propagating 
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undisturbed when the cloak is active. The comparison for a soft cyhnder is shown in Figure 




Figure 15: Absolute value of total pressure field when cloaking devices are inactive (left) and 
active (right) for scattering from the hard cylinder. Calculations are performed for a hard 
cylinder with M = 5 active sources, angle of incidence = 17°, and wave number ka = 5. 




Figure 16: Real value of total pressure field when cloaking devices are inactive (left) and 
active (right) for scattering from a hard cylinder. 



5 Discussion 

5.1 The case of many sources: M ^ 1 

The numerical results of §4.21 indicate better convergence properties for fewer multipoles 
if more active sources are used. It is therefore of interest to consider the limit in which 
many sources are available: the large M limit. Staying with the configuration of §4.11 the 
appropriate limit to consider is (a^ =)a ^ bn/M so that fca <C 1 and the small argument 
approximation can be used for the Bessel functions Jn{kam) (note kb is not necessarily small). 
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Figure 17: Absolute value of total pressure field with cloaking devices are inactive (left) and 
active (right) for scattering from a soft cylinder: k = 5^ ip = 17° . 

This implies that to leading order in ka the coefficients in (l7bl) reduce to 

r(-2)[/+'(xj, / = o, 

b^,in = ^ka X I le-'''- U+{^J, I = ±1, (34) 
[o, l^O,±l. 

The identity J^_i(x) — J^+i(x) = 2J'^{x) has been used to simplify the / = term in ( l3ll) . 
The source field follows from eq. f l2cl) and the identity l^Ji(x) = —V-^{^m) as 

oo M 

n=—oo m=l (35) 

[f/+(x^)//J'\|x - x^l) COS ( arg(x - x^) - 9^) - t/+'(x^)//i'H|x - x^|)' . 

The field of the active sources is therefore composed of monopoles and dipoles only, with 
no contribution from higher multipoles. This agrees with what one might expect from the 
continuous limit of M ^ oc, i.e. a closed contour of monopoles and dipoles, but here it is 
obtained from the discrete solution. In fact, eq. ( l35l) is 

M 

^d= ^CL^ K(Xm)5nVb(x - X^) - l/o(x - X^)^^^/^ (x^)] , (36) 
m=l 

which can be seen to be the discretized version of the fundamental integral identity eq. 
f l20l ) with the correspondence JdS ^ 2a While eq. f l34l ) is thus the natural first 

approximation for Ud based on the integral equation f l20l ). it should be realized that it was 
obtained here as a first order approximation of the exact expression f |7bl ). The latter therefore 
provides the basis for a multipole expansion of the exact source field obtained by including 
higher powers of ka than considered in eq. f l34l ). This possibility goes beyond our present 
interests but will be examined in a separate study dealing with approximations to the exact 
results of Theorem [2l 
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5.2 Zero scattering matrix 

The exact source field Ud{'x.) of eq. f l2cl) exhibits some interesting features. This field is, 
by design, equal to the negative of the incident field in the cloaking region C, and it also 
vanishes identically outside the concave region R defined in (Ej) The non-radiating property 
of Ud is as important as the fact that it exactly cancels the incident wave in C. Let us 
examine this more closely. Define the infinite matrix S with elements Spq such that 



M 



pq 



EE 



bm.lqUp 



m=l / 



or, usmg 



(37) 



M 



ka. 



Ln=—oo 



l + n 



Am) 



Am) • 



(38) 



The matrix S is, formally at least, like a scattering matrix. For instance, by inspection, 
S is hermitian [Spq = 5'*^). However, by design and based on Theorem [H S = 0, and 
as such it could be called a zero-scattering matrix. Alternatively, it can be viewed as a 
formula for generating non-radiating fields. This has relevance to the inverse source problem 
pH] . It is known that solutions to the inverse source problem are non-unique [I2], although 
some uniqueness results are available for restricted forms of sources, e.g. "minimum energy 
sources" [H]. The solution of the active cloaking problem as developed here has generated 
a new family of non-radiating sources, with the property that they cancel a given incident 
field over a finite region. 



6 Conclusions 

By definition, an active source cloaking strategy requires solution of an inverse problem: 
find the active source amplitudes associated with a given incident field in order to exactly 
cancel the latter in some finite region. The results given in Theorem [2] provide closed- form 
solutions for the inverse problem for an arbitrary time harmonic incident wave field. These 
new expressions require only the expansion of the incident field into entire cylindrical waves 
and can be evaluated to any degree of accuracy by increasing the truncation parameter 
associated with the number of modes of the active source. Simultaneously the fact that the 
active source field has been shown to vanish identically outside the region R defined in (|8j) 
means that the active field is non-radiating. This latter property is just as important as its 
ability to nullify the incident wave in the region C. 

The necessary and sufficient conditions on the active source coeflScients, given in Theorem 
[H provide a means to quantify the error in active cloaking when the number of modes is finite. 
These errors have been analyzed here in some specific scenarios. It has been shown that 
the error in the far-field amplitude decreases as N increases, M increases and k decreases. 
In particular there is a great sensitivity to the increase in A; relatively small errors can be 
attained in the far-field amplitudes for moderate A, say A ^ 10. On the other hand for small 
errors in the near-field amplitudes, relatively large values of A are required. Furthermore, 
there is a striking reduction in error when moving from the case of M = 3 to M = 4 
motivating the latter as a preference. In contrast to the far-field case, errors decrease for 
increasing k. 



Active exterior cloaking 



July 23, 2012 



20 



Numerical results were given which illustrate the cloaking effect in various instances, in- 
cluding the presence of a sound-soft and sound-hard circular cylinder. In the appropriate 
limits, perfect theoretical active cloaking is achieved. The availability of closed- form ac- 
tive source amplitudes opens the door for possible studies on practical realization of active 
cloaking devices. 

The case of many sources, where the active field degenerates to one involving a sum of 
monopole and dipole sources is worthy of further, separate study relating to the multipole 
expansion associated with the active field. Finally, the non-radiating nature of the active 
source field is especially noteworthy. The associated scattering matrix, defined in f l38l ) (which 
is zero by design) is therefore associated with a new family of non-radiating source solutions 
which would appear to be useful in the so-called inverse source problem. 
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